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Foreword 

Mathieu functions have many applications, especial- 
ly in mathematics and physics: 

• separation of variables in elliptical coordinates 

• separation of variables in parabolic coordinates 

• vibrations in a stretched elliptical ring membrane 

• stability of a pendulum with periodically varying 
length 

• amplitude distortion in moving-coil loudspeakers 

• relativistic oscillators 

These applications are mentioned in the chapter on 
Mathieu functions, written by Gerhard Wolf for the 
NIST Digital Library of Mathematical Functions. 
Professor Wolf is coauthor (with J. Meixner and F. W. 
Schafke) of Mathieu Functions and Spheroidal 
Functions and their Mathematical Foundations, Lecture 
Notes in Mathematics 837, Springer- Verlag, 1980. 

The DLMF is scheduled to begin service in 2009 
from a NIST Web site. A hardcover book will be pub- 
lished also. These resources will provide a comprehen- 
sive guide to the higher mathematical functions for use 
by experienced scientific professionals. 



The DLMF is modeled after the enormously success- 
ful but increasingly out-of-date NBS Handbook of 
Mathematical Functions, National Bureau of Standards 
Applied Mathematics Series 55, M. Abramowitz and I. 
A. Stegun (editors), 1964. The NBS handbook has sold 
more than 700,000 copies and is frequently cited in sci- 
entific journal articles. The need for a modern reference 
is being filled by NIST editors and staff, aided by the 
scientific content provided by approximately 50 exter- 
nal authors and validators 

In addition to the main purpose of the DLMF, which 
is to provide a comprehensive and authoritative 
research tool, the project also seeks to guide further 
research in special functions. The paper that follows is 
an example. It provides the proofs of results that Wolf 
presents for the first time in §28.4(vii) of the DLMF. 

Daniel W. Lozier 

DLMF General Editor 

NIST Mathematical and Computational Sciences 

Division 
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1. Introduction and Definitions 

In this paper we aim to give asymptotic formulae for 
Fourier coefficients of the periodic solutions of 
Mathieu's equation 



w + (a-2qcos2z)w = 0. 



(i) 



Equation (1) possesses the fundamental pair of solu- 
tions Wi(z; a, q), w n (z; a, q) called basic solutions (see 
Ref. [l])with 



WjiOi^q) w E (0;a,q) 
Wj(0;a,q) w^(0;a,q) 



1 
1 



(2) 



Furthermore, we obtain eigenvalues and eigenfunctions 
of(l) for #2 = 0, 1,2,.... 

Table 1 gives their notations and properties. "Period 
if means that the eigenfunction has the property 
w(z + 71) = w(z), whereas "Antiperiod 7i" means that 
w(z + n) = -w(z). "Even parity" means w(-z) = w(z) and 
"Odd parity" means w(-z) = -w(z). 

The Fourier coefficients satisfy the recurrence rela- 
tions 

H ~ a A = °> a = <hn(<?)> Am = 4m(4)> 

(a-4)A 2 ~q(2A +A 4 )=0, 

(a-4m 2 )A 2m -q(A 2m _ 2 +4 m+2 ) =0, m =2,3,4,..., 

(3) 

(a-l-q)A 1 -qA 3 =0, a = a 2n+1 (q), A, m+l =4^M\ 

(a-(2m+l) 2 )A 2m+1 -q(A 2m _ 1 +A 2m+3 ) =0, m=l,2,3,..., 

(4) 



(a-4)B 2 -qB 4 =0, a = b 2n+2 (q\ B 2m , 2 = %£{q) 
(a-4m 2 )B 2m -q(B 2m _ 2 +B 2m+2 ) = 0, m=2,3,4,.... 



(6) 



In §2 we examine the asymptotic behavior of the coef- 
ficients At{q\ A™(q), B\^ y {q\ and B%*(q) for fixed 

77, q (^ 0), as m — » ©°. 



2. Asymptotic Forms 

The following result will be proved: 
Proposition 1. For fixed n and q^Q, as m — > °° 

4^ (?) _H)" f*Y wa+oa/«)) 



4\q) (m!) 2 U I w^\n;a 2 M\q) 



(V) 



4TM _ (-1)- +1 (q 



4" +1 (q) (0 m+1 ) 2 l 4 J <(iK>"2„ +l (q\q) 



2(l + 0(l/w)) 



-. (8) 



B? m + Mq)_ (-D" (q 



*r(q) ((*)..♦,) (4 ) w^n;b 2n ^q\q) 



2(1 + 0(1 lm)) 



; (9) 



BIT (?) _ H) m ( q T «*0 + ofl / «)) 



-• (10) 



B^fo) (th!) 2 ^4J H^^fo),*) 

Proof of (7). We consider for ^ ^ the three term- 
recurrence relation 



(2m) 2 



Q 



\6m 2 (m-iy 



' Z m-1> w=2,3,. 



(11) 



(a-l + gO-8, -^ 3 -0, a-fe^g), 5 2m+1 -Ct!(^)> For the two independent solutions u m v m of (11) with 

(a-(2w+l) 2 )5 2m+1 -q(B 2m _, +B 2m+3 ) =0, m=l,2,3,.~, u l = a,u 1 = a (1 -\a) +\q 2 and v, = 1, v 2 = (1 - ±a), it 



(5) 



Table 1. Eigenvalues and eigenfunctions 



Eigenvalues 


Eigenfunctions 


Periodicity 


Parity 


Fourier series 


a = a 2n (q) 


ce 2fl (z, q) 


Period 7T 


Even 


XI =0 4:(<7)cos(2/n)z 


a = «2iH-l(?) 


ce 2/J+1 (z, q) 


Antiperiod tt 


Even 


II = „4: + + 1 1 (9)cos(2m+l)z 


a = b 2n+1 (q) 


se^ife #) 


Antiperiod 7T 


Odd 


XI.o^ , ifa) sin ( 2m+1 )z 


a = b 2n+2 (q) 


se^O, #) 


Period n 


Odd 


H = o 5 2:: 2 2(9)sin(2/n+2)z 
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follows from Satze (Theorems) 1 and 3 of F. W. 
Schafke [2] that 



\imu m =-(2/7r)Wj((7r/2);a,q), 

m— >» 

limv m =(2/7r)w JI ((7r/2);a,q). 

m— >» 

Transformation of (11) with 



(12) 



C 2m =(-l) m - 1 —((m-l)]) 2 z M , m=l,2,3,..., (13) 



yields 

^( C 2( W+ l)+ C 2( m -l)) = («- 4 ^ 2 )Qm ? W=2,3,.... (14) 

We note the following special solutions of (14): 

U 2m with z m = u m , V 2m with z m = v m . (15) 

Then 



U 2 =a, U 4 =-(4/q)(a(l-±a)+±q 2 ), 
V 2 =l, V 4 =-(4/q)(l-±a). 



(16) 



Furthermore, we have 

U 2im+1) V 2m - V 2(m+1) U 2m =U 4 V 2 -V 4 U 2 = -2q. (17) 

If, now, a = a 2n {q), then w\{tzI2\ a^iq), q) = and u m — > 
as m —> °°. 

Set C/ = #. Then for m = 2, 3, 4, ..., 



aU -qU 2 =0, (a-4)U 2 -q(2U +U A )=0, 

q{U 2im+l)+ U 2{m _ X) )={a-Am 2 )U 2m . (18) 



Comparison with^(^) of (3) shows that 

Vfo) 2m ' 



(19) 



Thus U 2m is the minimal solution of (14), and by the 
substitution 



r = l-*l -Up, 



1 4 J (w!) 2 '' 



(20) 



we find that 



j 02. (g) 

(2m) 2 



IP- -P— i =77 2, .^ P-ti' ( 21 ) 
16 w (w+1) 



P„ " P„-i = 0(1 / m 2 ), and p m = k + O (Vm). (22) 

The constant k is determined with the aid of (17): 
it V -V I J 

U 2(m+l) y 2m y 2(m+l) U 2m 

\ m+1 1 J_m-1 

1) «^ljiF' , -« < - 1) "V << " , - 1)!)2v - 



*T i 



4 1 (w!) 2 



-P m+ iv m -v m+1 p m =-2q, (23) 



4 1 (m(m-l)) 



and for m — > °° 

(2t/wK(iff;a l ,( 9 ), 9 )=2* (24) 

Together with (19) we obtain (7). 

Proof of (10), In the same way it follows, if a = 
K+iiq), then w n (j7t; a, q) = 0, v m -> 0, and 

(*-4)*W 4 =0, q(V 2(m+1) +r 2(m ., } ) = (a-4m 2 )^ m , 
m = 2,3,4,.... (25) 

Comparison with B%£ 2 (q) of (6) shows that 

j2n+2 / 



^rw 



2n+2/ \ 2/m- 



Thus K 2m is the minimal solution of (14), and for 



q\ l 



(26) 



(27) 



4) (m}) 2 
we find again p m - p m _ 1 = 0{\im 2 ) and 

p m =k + 0{\lm). (28) 

The constant k can be computed via (17): 



tj V -V U 

,gt 1 

= " m+ iP m -| t I 7-7— tttP». + i"». 



4 1 (m(m + \)Y 



= -2q. 
Letting m —> ©°, we obtain 



(29) 
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qn 



Wi(hn\KA<i)><l) 



(30) 



Together with (26) we obtain the formula (10). 

Proof of (8) and (9). We start with the recurrence 
relations 



(2m + 1) 2 



\6{m+\)\m-\y 



-^ w _i, m =1,2,3,.... 



(31) 



The two independent solutions u m and v m of (3 1 ) with 
u = 1, u x = q - a + 1 and v = 1, v 1 = -q - a + 1, respec- 
tively, satisfy 



limw m =w 1 (K/2;a,q), 
limv m = w' a (nl2;a,q); 



(32) 



see Satze (Theorems) 1 and 3 of F. W. Schafke [2]). 
Then transformation of (31) with 1 

C 2m+1 =(-4/g) m g((^) m ) 2 z m , m =0,1,2,3,..., (33) 

yields 

?(C ta+3 + C 2m - l ) = («-0 +2w) 2 )C 2m+1 , m =1,2,3,... . 

(34) 

We note the special solutions 

U 2m+l with z m = m m , K 2m+1 with z m = v m , (35) 

and obtain 

^=4, U 3 =a-l-q, V l =q 9 V 3 =a-l + q (36) 

Furthermore, we observe that 

U 2m+1 V 2m _, - V 2m+l U 2m _, = £/ 3 ff - ^ = -2*>. (37) 

If, now, a = a 2 «+i(#)> then w I (^/2; d^-iteX #) = and 
u m -> for m -> <*>. 

Comparison with ^+\(<7) of (4) shows that 



(-j) m means the Pochhammer symbol: (a) n = a(a + \){a + 2) ... (a 



% 2 " +1 (<?) 2m+1 ' 



(38) 



Thus £4«+i is the minimal solution of (34), and with 



u 2m+1 = 



q\ 1 



(0- + i) 



(39) 



we find that 



j_ a 2„ + i(<?) L . _ 



(2m + \y 



Pm-Pm-I 



(2m+l) 2 (2m+3) 2 Pm+1 ' 
(40) 



Pm-p m -^0(\im 2 \ p m =k + 0(l/m). (41) 

The constant k is determined with the aid of (37): 

U^V^-V^U^ =(Kl/m 2 )-qv mPm _ l =-2<f. 

(42) 



Letting m — > <*>, we obtain 

2<7 



* = 



Wn(^/2;a 2fl+1 (^) 9 ^)' 



(43) 



Together with (38) we find (8). 

The formula (9) is obtained by the transformation 
q -> -<?, a^HO = *2«+i(?), and A^(-q) = 



3. Improvement of the Rate of 
Convergence 

If a^ ^ (2m) 2 , then for m = 0, 1 , 2, ... we transform 
(21) via 



p.=wnh 



(2r) 2 



(44) 



and obtain (as in Ref. [2]) t, m - 1, x = Q {\lm\ £,„-£,= 
0(]/m 3 ), where £ = iff sin (^n^a 2n (q) )fc. Finally we 
have 

^(g) = (-i)"fgT ? , 

^fo) (tk!) 2 ^ 

2sin(^7^))(l + 0(l/m 3 )) 



*-> (45) 
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In the other cases, for example (40), similar results 
can be obtained. 

Error bounds are possible with the aid of Refs. [2] 
and [3]. 

Remark: Further improvements of the rate of con- 
vergence are possible by application of the results of 
Ref [3J. 
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